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Abstract 

A  quantum  lattice  algorithm  is  developed  to  examine  the  effect  of  an  externa,  potential  well  on  exactly  integrable  vector 
Manakov  solitons.  It  is  found  that  the  exact  solutions  to  the  coupled  nonlinear  Schrodinger  equations  act  like  quasi- 
solitons  in  weak  potentials,  leading  to  mode-locking,  trapping  and  untrapping.  Stronger  potential  wells  will  lead  to  the 
emission  of  radiation  modes  from  the  quasi-soliton  initial  conditions.  If  the  external  potential  is  applied  to  that  particular 
mode  polarization,  then  the  radiation  will  be  trapped  within  the  potential  well  The  algorithm  developed  leads  to  a  finite 
difference  scheme  that  is  unconditionally  stable.  The  Manakov  system  in  an  external  potential  is  very  closely  related  to  the 
Gross-Pitaevskii  equation  for  the  ground  state  wave  functions  of  a  coupled  BEC  state  at  7’  =  OK. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  optical  propagation  down  a  birefringent  fiber,  under  appropriate  conditions  [1]  the  two  orthogonal 
polarization  amplitudes  satisfy  the  coupled  nonlinear  Schrodinger  equations  (NLS) 

iSrGi  +  9**01  +  2(10,  |2  +  5|02|2)0,  =  0  , 

'9(02  +  0**02  +  2(1  Q2 12  +  B\Qx\2)Q2  =  0 ,  (1) 

where  B  is  the  cross-phase  modulation  coefficient  (XPM):  B  =  (2  +  2  sin2  0)/(2  +  cos2  8)  and  6  is  the 
birefringence  angle  of  the  fiber.  For  B  =  1,  Eq.  (1)  is  exactly  integrable  by  the  inverse  scattering  as  well  as  the 
Hirota  method  to  yield  the  so-called  Manakov  vector  solitons.  Vector  solitons  have  a  richer  behavior  than 
their  scalar  counterparts  that  remain  invariant  in  amplitude  and  speed  under  soliton  -soliton  collision.  The 
only  signature  of  a  scalar  soliton  collision  is  an  induced  phase  shift.  However,  in  vector  2-soliton  interactions, 
depending  on  the  initial  amplitudes,  there  exists  the  possibility  of  inelastic  collisions  [2], 
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Polarization  #1  -  MANAKOV  SOLITONS  Polarization  #2  -  Manakov  Solitons 


0.000  0.010  0.020  0.000  0.025 


Fig.  1.  Long  time  evolution  of  the  Manakov  soliton  amplitudes  12,1  and  |g2|  under  periodic  boundary  conditions.  Time  increases 
vertically  down.  The  total  intensity  in  each  mode  is  numerically  preserved  by  our  algorithm  to  better  than  1  part  in  107.  All  soliton  speeds 
are  equal  and  are  invariant  under  collisions.  Soliton  amplitudes  are  invariant  in-between  collisions  and  can  be  changed  by  a  collision, 
subject  to  conservation  of  the  total  intensity  in  each  mode. 


If  one  employs  periodic  boundary  conditions  (which  permits  many  Q,  -  Q2  soliton  collisions),  there  is  on 
the  next  collision  the  re-appearance  of  the  soliton  destroyed  in  the  first  collision  as  can  be  seen  in  the  long  time 
evolution  of  the  Manakov  solitons  Fig.  1 . 

The  4  soliton  speeds  are  time-invariant,  with  only  their  amplitudes  being  affected  by  the  collisions  subject  to 
constant  total  intensity  in  each  mode.  For  the  parameters  chosen  the  soliton  collisions  always  occur  at  the  end 
point(s)  and  at  the  mid-point  (on  ignoring  the  very  small  collision-induced  phase  shifts).  Furthermore,  the 
corresponding  right-traveling  solitons  from  each  polarization  are  mode-locked  and  similarly  for  the 
corresponding  left-traveling  solitons:  i.e.,  the  center  of  mass  of  each  right-traveling  soliton  in  the  two 
polarizations  are  coincident  (as  are  the  center  of  masses  for  the  left-traveling  solitons). 

Here  we  shall  consider  the  effect  of  slowly  (spatially)  varying  external  potentials  on  this  vector  Manakov 
system 

i0/2i  +0~e,  +2(ie,|2  +  |<22I2)Gj  +  V?\x)Q ,  =  0  , 

10102  +  0XX02  +  2(1  g2 12  +  16,  |2)02  +  V?\x)Q2  =  0  .  (2) 

For  nonzero  external  potentials,  the  Manakov  soliton  solutions  of  Eq.  (1)  can  at  best  be  quasi-solitons  since 
Eq.  (2)  is  now  nonintegrable.  It  is  precisely  this  that  we  shall  investigate  here,  keeping  in  mind  that  Eq.  (2) 
bears  a  striking  resemblance  to  the  Gross-Pitaevskii  [1]  equation  that  describes  the  ground  state  wave  function 
for  the  Bose-Einstein  condensate  (BEC)  state  at  T  =  0  K  in  an  external  confining  potential. 


2.  Vector  Manakov  equations  in  external  potentials 

A  quantum  lattice  representation  for  the  vector  Manakov  system  Eq.  (1)  has  been  presented  in  detail  in  Ref. 
[3]  so  we  shall  only  briefly  describe  the  algorithm  here.  After  discretizing  the  spatial  domain- one  introduces  2 
qubits/node  for  each  polarization  mode.  Since  the  diffusive  part  of  Eq.  (1)  does  not  couple  the  polarizations 
Qx  and  Q2,  this  can  be  recovered  from  the  quantum  lattice  algorithm  by  simply  collisionally  entangling  just 
the  2  qubits  of  the  same  polarization  mode  at  each  spatial  node.  The  required  unitary  collide-stream 
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sequence  is  [3] 

(  0\U  +  At) 

[  Qz(.t  +  At) 


[SlCS:CST CS:C.SjCS,  csjcs,  C]Q,(r) 
[Sj  CS4  CSj  CS4  C.s]  CS3  CS]  CS}C]Q2(t) 


where  the  global  unitary  collision  operator  C  is  a  tensor  product  of  the  local  VSWAP  gate  that  entangles  the  2 
qubits  of  each  polarization  separately,  while  the  streaming  operator  on  qubit  “i”  spreads  this  entanglement 
throughout  the  lattice.  5,  streams  qubit  “i”  to  the  right  while  Sj  streams  qubit  “i”  to  the  left.  The 
symmetrizing  of  the  algorithm  on  qubits  “2”  and  qubits  “3”— “4”  will  result  in  a  higher-order  accurate 
finite  difference  scheme.  In  the  continuum  limit  under  diffusive  ordering,  Ax2  =  0(s2)  —  At,  Eq.  (3)  yields  the 
uncoupled  free-Schrodinger  equations  for  the  polarization  modes.  (Ax  is  the  lattice  cellsize)  [4].  The  nonlinear 
coupling  of  the  polarizations  is  readily  obtained  by  introducing  appropriate  phase  modulations  into  the 
polarization  wave  functions  0,  and  Q2.  This  can  be  understood  from  path  integral  formulations  of  quantum 
mechanics  or  from  multiscattering  theory  of  classical  electromagnetic  wave  propagation  [5]  in  which  phase 
modulations  are  introduced  and  then  converted  into  amplitude  modulations  by  the  diffusive  part  of  the 
parabolic  Schrodinger  equation.  Thus,  with  the  phase  modulation 


applied  following  the  unitary  collide-stream  operator  sequence,  Eq.  (3),  one  obtains  in  the  continuum  limit 

+3 **0,  +  +  0(e2)  =  0  , 

i3/£?2  +3**02  +  V2(\Q\\,  IQ2O82  +  ^(s2)  =  0  ,  (5) 

where  the  potentials  V,-  (i  —  1 , 2)  are  arbitrary  functions  of  the  polarization  wave  functions  and/or  space.  A 
particular  choice  of  these  potentials  thus  yields  the  required  vector  Manakov  equations  inn  external  potentials, 
Eq.  (2). 

3.  Manakov  quasi-solitons  in  an  external  potential  well 

Taking  as  initial  conditions  the  Manakov  solitons  of  Fig.  1,  the  role  of  an  external  potential  well  of  depth 
t/o  on  the  subsequent  time  evolution  is  considered,  with 


V extW  =  U0  tanh 


fx-Lx/2\ 

v  Lx/ 5  r 


0<x<E* 


Vext  is  symmetric  about  L*/2  where  it  is  zero,  and  then  gently  asymptotes  to  t/o  as  x  ->  0  or  Lx. 

We  first  consider  the  effect  of  a  very  weak  potential  well  on  the  polarization  mode  0,  only:  i.e.,  Eq.  (2)  with 
Kft(x)  =  Fext(x),  and  k^xt(x)  =  0.  Here,  the  potential  well  t/o  is  an  order  of  magnitude  smaller  than  the 
amplitude  of  the  smallest  initial  soliton  in  Fig.  1.  The  external  potential  is  sufficiently  weak  that  the  initial 
conditions  yield  quasi-solitons  that  retain  their  identity  under  collisions  and  there  is  mode-locking  of  the 
quasi-solitons.  However,  the  external  potential  creates  enough  perturbation  that  both  quasi-solitons  of  mode 
0,  now  still  exist  after  the  first  collision.  Since  the  quasi-solitons  are  mode-locked,  one  need  only  present  the 
long  time  evolution  for  one  of  the  polarization  modes.  In  Fig.  2,  we  plot  the  evolution  of  the  quasi-solitons  of 
mode  02.  The  initially  right  traveling  quasi-soliton  undergoes  an  amplitude  modulation  after  the  1st  collisipn 
and  then  becomes  trapped  in  the  weak  potential  well  (/  =  50  x  4.5  k).  It  then  moves  to  the  left,  parallel  with 
the  other  soliton  (which  is  always  untrapped)  and  reflects  off  the  right  potential  wall  {t  =  105  x4.5k). 
However,  after  the  next  collision  (around  t  =  150  x  4.5  k)  this  quasi-soliton  becomes  un trapped  and  stays 
untrapped  for  all  times. 

The  same  very  weak  external  potential  is  now  applied  to  mode  2,  with  no  external  potential  on  mode  1: 
V\n(x)  =  0,  and  F^xt(x)  =  Fext(x).  Because  the  initial  state  for  mode  2  is  very  different  from  that  for  mode  1 , 
the  effect  of  the  external  potential  on  mode  2  will  be  different  from  that  on  mode  1.  The  run  time  is  now 
extended  to  3.6 M  iterations  and  is  plotted  in  Fig.  3  in  two  parts:  0<f<  1.8 M,  and  1.8M<?<3.6M.  Again, 
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Poiorizotion  #2 


X  *  16 


0.000  0.026 
VI  =  0.001 ,  <knem=  2900 

Fig.  2.  The  effect  of  a  very  weak  external  potential  well  on  mode  1 .  The  quasi-solitons  are  mode-locked  so  only  the  time  evolution  of  mode 
2  is  presented.  After  the  1st  collision  the  right  traveling  quasi-soliton  is  trapped  by  the  potential  well,  but  it  escapes  from  the  well  after 
basically  one  period.  The  initial  left-traveling  quasi-soliton  is  always  untrapped. 


Polarization  #2  Polarization  #2  s  TIME  *  1.8  M  ~>  3.6  M 


0.000  0.025  0.0100.0200.030 

V2  =  0.001 ,  denow=  2500  V2  =  0.001 ,  denom=  2500 


Fig.  3.  The  effect  of  the  same  very  weak  external  potential  well  as  in  Fig.  1,  but  now  only  on  mode  2.  The  quasi-solitons  are  mode-locked 
so  only  the  time  evolution  of  mode  2  is  presented,  but  this  time  for  twice  the  number  of  iterations  plotted  in  Fig.  2.  It  is  now  the  initially 
left-moving  quasi-soliton  that  is  trapped  and  remains  trapped  by  the  external  potential  for  all  time.  After  4  collisions,  the  initially  right 
moving  soliton  becomes  trapped  (around  t—  340  x  4.5  k)  for  one  period,  and  then  becomes  untrapped  and  then  eventually  retrapped 
(i=  1.8  M  +  300x4.5  k). 


the  quasi-solitons  are  mode-locked  and  so  results  are  presented  only  for  mode  2.  The  left-moving  quasi-soliton 
at  t  =  0  remains  trapped  by  this  very  weak  potential  well  for  all  time,  0<f <3.6M.  On  the  other  hand,  the 
t  =  0  right-moving  quasi-soliton  undergoes  several  untrapped-trapped  sequences. 
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The  strength  of  the  potential  well  is  now  increased  by  a  factor  of  2.5.  Again  the  resulting  quasi-soliton 
motion  is  mode-locked  and  we  present  the  results  on  mode  polarization  1.  For  early  times,  0<t <265  x  4.5  k, 
the  potential  well  is  sufficiently  deep  that  both  initial  quasi-solitons  are  trapped.  However,  for  later  times  there 
is  an  intricate  interwined  sequence  of  trapped-untrapped  that  occurs  among  the  quasi-solitons,  as  seen  in 
Fig.  4.  If  the  potential  well  is  now  applied  to  the  polarization  mode  Q2  rather  than  to  Q{,  the  two 


Polarization  #1 


0.000  0.010  0.020 

VI  =  0.0025  ,  <fcnom=  2500 


Polarization  #1 :  t  =  1.8  M  ==>  3.6  M 


0.000  0,010  0.020 

VI  =  0.0025 ,  dcnoiip  2500 


Fig.  4.  The  effect  of  increasing  the  potential  well  by  a  factor  of  2.5.  The  quasi-solitons  are  still  mode-locked  for  polarizations  1  and  2.  With 
the  same  initial  conditions  for  all  runs  reported  in  this  paper,  we  now  find  that  the  quasi-solitons  are  initially  trapped  and  after  more 
collisions  an  interweaving  alternating  web  between  the  two  quasi-solitons  of  trapped-untrapped. 


Polarization  #1 


0,000  0.010  0.020 

V2  =  0.0025,  <tenom=  2500 


Fig.  5.  The  same  external  potential  as  in  Fig.  4,  but  now  applied  to  the  polarization  mode  Q2.  Again  there  is  mode-locking  between  the 
two  polarization  modes  2i  and  Q2.  The  evolution  of  the  quasi-solitons  for  mode  £>,  are  plotted  up  to  time  t=  1 .8  M,  and  both  quasi- 
solitons  are  trapped  within  the  potential  well.  The  difference  between  Figs.  4  and  5  arise  from  the  application  of  the  external  potential  to 
different  polarization  mode  initial  quasi-soliton  amplitudes. 
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quasi-solitons  are  trapped  for  all  times,  shown  in  Fig.  5.  The  external  potential  is  now  increased  by  a  factor  of 
5  over  that  in  Fig.  3.  It  is  applied  to  polarization  mode  Q2  only.  In  Fig.  6  we  plot  the  evolution  of  the  initial 
quasi-solitons  for  both  mode  (with  no  external  potential)  and  mode  Q2.  For  times  t<  1 70  x  4.5  k,  there  is 
mode-locking  between  the  <2i  and  Q2  polarizations.  Shortly  after  this  one  starts  to  see  radiation  modes 
between  spewed  out  in  the  Q{  polarization  and  eventually  these  modes  are  unconfined  (Fig.  7  for 
t  =  300  x  4.5  k).  However,  for  polarization  Q2  radiation  modes  are  emitted  later  in  time  but  they  are  confined 
by  the  potential  well  applied  to  Q2  (Fig.  7b  for  t  —  300  x  4.5  k).  For  r >  350  x  4.5  k,  quite  strong  radiation 
modes  are  emitted  in  the  Q2  polarization,  but  they  are  confined  by  the  potential  well.  As  the  external  potential 
well  is  increased,  so  will  the  potential  gradients  until  eventually  the  t  =  0  quasi-solitons  are  completely 
destroyed  and  only  radiation  modes  exist.  If  the  external  potential  is  applied  to  polarization  Qt  then  the 


Polarization#!  Polarization  #2 


0.000  0.010  0.020  0.00  0.01  0.02  0.03 
V2s  0.005 ,  denom=  2500  V2  =  0.005 ,  denom=  2500 


Fig.  6.  The  mode  locking  of  Q,  with  g2  is  broken  for  t  >  170  x  4.5  k,  with  radiation  emission  in  both  Q,  and  Q2.  For  g,  the  radiation  is 
untrapped,  while  trapped  for  g2. 


polarization  #1 


X  *  16  X*1  6 

Fig.  7.  A  snapshot  at  time  t  =  300  x  4.5k  of  |gj|  and  |g2|  with  the  external  potential  only  applied  to  polarization  g2.  The  quasi-soliton 
structures  still  exist  and  are  mode-locked  for  the  polarizations  (located  around  x  =  150  x  16  and  x  =  375  x  16).  However,  radiation 
modes  are  also  emitted  in  both  the  g[  and  g2  solutions.  In  the  case  of  g1  the  radiation  occurs  throughout  the  entire  lattice  while  for  g2  the 
radiation  is  spatial  confined  by  the  potential  well  (here  to  the  region  225  x  16<x<550  x  16). 
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radiation  modes  generated  with  the  collapse  of  the  initial  quasi-solitons  will  be  trapped  by  the  external 
potential. 

4.  Conclusion 

A  quantum  lattice  representation  for  vector  solitons  in  an  external  potential  has  been  developed.  The 
algorithm  employs  a  sequence  of  unitary  collision  operators  (VSWAP  gate)  that  locally  entangle  the  qubits  at 
each  node  and  unitary  streaming  operators  that  rapidly  spreads  this  entanglement  through  the  lattice.  It  will 
be  difficult  to  implement  our  quantum  algorithm  experimentally.  First,  the  quantum  phase  coherence  must  be 
maintained  throughout  the  complete  sequence  of  collide-stream,  Eq.  (3).  This  must  then  be  followed  by 
quantum  measurement  from'which  we  need  to  extract  the  absolute  values  of  the  amplitudes  and  not  simply  the 
probabilities  which  are  naturally  extracted  from  the  projective  von  Neumann  quantum  measurement  over  a 
large  ensemble.  This  is  how  the  nonlinear  interaction  potentials  are  determined,  Eq.  (4).  Our  quantum  lattice 
algorithm  gives  rise  to  a  finite  difference  scheme  that  is  explicit  and  unconditionally  stable.  Moreover,  the 
quantum  algorithm  is  very  efficiently  parallelized  on  a  classical  computer. 
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